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ABSTRACT 



\ We report simultaneous multicolor near-infrared (NIR) observations of the 

supernova associated with x-ray Flash 060218 during the first 16 days after the 
high energy event. We find that the light curve rises and peaks relatively fast 



compared to other SN Ic, with the characteristic broad NIR peak seen in all 
three bands. We find that the rise profile before the peak is largely independent 
of NIR wavelength, each band appearing to transition into a plateau phase around 
day 10-13. Since the light curve is in the plateau phase when our observations 
end at day 16, we can only place limits on the peak absolute magnitudes, but 
we estimate that SN 2006aj is one of the lowest NIR luminosity XRF/GRB 
associated SNe observed to date. The broad peaks observed in the JHK S bands 
point to a large increase in the NIR contribution of the total flux output from 
days 10-16. This evolution can be seen in the broad color and SED diagrams 
constructed using UBVRIJHK S monochromatic flux measurements for the first 
16 days of the event. Ultimately, a 10-day rise time would make SN 2006aj an 
extremely fast rise SN Ic event, faster than SN 1998bw and SN 2003dh, which 
combined with its underluminous nature, indicates a lower amount of 56 Ni ejected 
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by the progenitor compared to other XRF/GRB-SNe. Furthermore, the lack of 
significant color change during the rise portion of the burst points to little or no 
spectral evolution over the first 10 days of activity in the NIR. We note that the 
onset of the plateau phase seen in the NIR roughly coincides with the reports of 
circularization of the SN ejecta indicating that the fast rise time and subsequent 
light curve behavior may be connected to the asphericity of the explosion. 

Subject headings: gamma-rays bursts — Infrared: Observations 



1. Introduction 



There is now a growing body of evidence to suggest that most long duration, soft- 
spectrum gamma-ray bursts (GRBs) are associated with the core collapse of massive stars. 
This connection had long been suggested on theoretical grounds (Colgate 1968, Woosley 
1993) as well as due to the similarities between GRB and SNe energetics. As of early 2006, 
this connection had been solidified by the detection of a total of three GRBs (GRB 980425/SN 
1998bw, GRB 031203/SN 20031w , and GRB 030329/SN 2003dh) with direct spectroscopic 
connections to supernova events. The early optical spectra of these supernovae components 
for which high-quality spectra exist show remarkable similarities, typically deficient of hy- 
drogen and helium lines, a property consistent with Type Ic supernova. Furthermore, very 
broad absorption lines of O I, Ca II, and Fe II indicate high expansion velocities and a 
kinetic energy (Ek) going much higher than the values inferred in typical SN Ic events. 
Despite great interest in understa nding the GRB-SN connections, only a few events hay e 
been observed well in the infrared (ILipkin et al.ll2004l ; ICobb et al.ll2006l ; iBloom et al.l 12004 ). 
Since estimates of E K and the amount of 56 Ni produced in the explosion depend on the total 
bolometric output, IR observations become particularly important for our understanding of 
the fundamental properties of these events. 

On 2006 February 18 at 03:34:3 UTC the Swift spacecraft detected and localized 
XRF 060218 (jCusumano et al.l l2006al ). The burst was immediately recognized as an un- 
usual event when compared to other bursts observed by Swift. It was far longer than any 
GRB or XRF previously detected by the spacecraft, with a duration of roughly 2000 
seconds. It also exhibited a soft gamma-ray spectrum ( Barbier et al .1120061 ) and an unusual 
afterglow that brightened for the first ten hours after the explosion before trans itioning to 
a more common power law decay (jCusumano et al.l l2006bl : ICampana et al.l 120061 ) . The first 
i ndications of an un derlying; SN were noted by spectra taken roughly 3 days after the event 
( IMasetti et al.ll2006l) followed by reports of a rebrightening of the optical counterpart shortly 
thereafter (jD'Avanzo et al.ll2006l ). Like previously observed XRF/GRB-SN the spectra of 
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SN 2006aj showed continuum and broad-line emission with significant amounts of Fe II and 
Si II A6355 but no H or He emission. The SN 2006aj spectra best matches those of SN 2002ap 
and SN 1997ef at similar epochs, both broad line SN Ic event s but neither of whic h have been 
conclusively shown to be associated with a XRF or GRB (IModjaz et al.ll2006l ). Emission 
line measurements by iMirabal et al.l (j2006al ) gave a redshift of z = 0.0331 making it one of 
the closest XRF/GRB ever observed, second only to GRB 980425. Such a close distance, 
~ 140 Mpc, meant that XRF 060218 was similar to two other nearby XRF/GRB-SN events, 
GRB 980425 and GRB 031203, in that they were all significantly underluminous in the total 
amount of emitted ga mma-ray, x-ray, and ra dio energies when compared to the majority 
of cosmological bursts (jSoderberg et al.l 120061 ) . The V band peak of SN 2006aj occurred 
at roughly T = 10.0 day s (9.7 days in the rest frame) with a peak apparent magnitud e 
of M v = 18.7 ± 0.2 mag JModiaz et alihood : ISollerman etaJbood : IMirabal et aDl2006ah . 
making it both the fastest evolving and least luminous of the XRF/GRB-SN observed to 
date. Pre-burst imaging of the SN 2006a j field by SDSS provided its host galaxy magni- 
tudes ( jCool et al.ll2006l ; iHicken et al.ll2006l ). whi ch given its distance reveal a low-luminosit y 



dwarf galaxy as the host, with low metallicity flModiaz et alibood: Iwiersema et al.l l2007h 



similiar to the majo rity of other GRB host galaxies (IFruchter et al. 



2006 



Stanek et al. 



2006 



Kewlev et al.l 120061 ) . 



PAIRITEL observations of XRF 060218/S N 2006aj start e d on Feb 2 0th, well before th e 
observed Fband and J band peaks reported by lModjaz et al.l (120061 ) and I Cobb et al.l (120061 ). 
and continued for the next 14 days. Here in Paper I we report on the light curve properties 
and host-subtracted absolute magnitudes of SN 2006aj in the J, H, and K s infrared bands 
during these epochs. We reserve a m ore detailed presenta tion of the broadband SED evolu- 
tion during this period for Paper II (IModjaz et al.l 120071 ) . We present our observation and 
data reduction techniques and results in §2 and discuss the implications of our observations 
to the further understanding of the XRF/GRB-SN phenomena in §3. 



2. Data & Analysis 



All of the infrared data presented in this paper were taken with the fully automated 
PAIRITEL project @ located on Mt. Hopkins in Arizona. The 1.3m telescope contains the 
three Near-Infrared Ca mera and Multi- Objec t Spectrometer (NICMOS3) arrays formerly 
of the 2MASS project (ISkrutskie et al.l 120061 ) to simultaneously read near infrared J, H, 
and K s (1.2, 1.6, and 2.2 /im respectively) band images. Each image consists of a 256 x 



http: / /www. pairitel.org 



-4- 



256 array with a pixel scale of 2 arcsec/pixel and an integration time of 7.8 s per image. 
These individual images are dithered in order to correct for bad pixels, and mosaics are 
created by drizzling the images, producing a Nyquist-sampled image with a pixel scale of 1 
arcsec/pixel. The telescopes entire operating software has been built on the Python program- 
ming lang uage and allows for autonomous operation which can be monitored and controlled 
remotely (IBloom et al.l 120051 ) making it excellent for the followup of transient events. 



Observations of XRF 060218 began on the morning of Feb 20th and continued until 
the object fell below our airmass limits some 16 days later. A total of > 9000 individual 
images were taken in the J, H, and K s wavebands to produce a total of 96 reduced mosaics. 
The individual images were dithered through a predetermined pattern for a single epoch of 
observations and then reduced through a custom pipeline. Bias and sky frames for each 
image were produced by median-combining several dithered exposures before and after the 
frame, whereas archival frames were used for th e bad-pixel masks and fl at fields. The final 
mosaics were produced using a drizzle technique (IFruchter fc Hooklll997l ) for each epoch and 
filter with an equivalent exposure time ranging between 2 and 30 minutes. The average 
exposure time of each mosaic generally increased as the object moved to higher air mass in 
subsequent nights. Of the 96 mosaics produced through the automated reduction pipeline, 
3 were unusable due to poor transmission on three nights, leaving ~ 31 usable mosaics per 
filter. 

A custom pipeline was used to perform photometry of the final J, H, and K x mosaics, 
using aperture photometry via the Sextractor package (IBertin &: Arnoutd Il996l ) and PSF 
photometry via the NSTAR routine in IDL to estimate the instrumental magnitudes of 
every object in each mosaic. These values were then compared to the original 2MASS 
catalog of the same field for zeropoint determination of each reduced frame. Finally, an IDL 
routine was employed to produce a light curve for every stellar object in the field for a final 
relative calibration. The resulting median Am variations in these stellar light curves allow 
for additional corrections to be made to the SN light curve that account for any errors in 
the zeropointing of the individual mosaics. 

Overall, we find that the aperture and PSF fitting routines generally yield equivalent 
values for J, H, and K s band mosaics. We ultimately chose to use aperture photometry with 
a floating aperture radius roughly equal to 1.5, 1.75, and 2.0 times the seeing (FWHM) for 
the J, H, and K s bands respectively. These factors were empirically determined to maximize 
the signal to noise in each band and result in aperture radii roughly ranging from 3 to 5 
pixels. 

To account for the host galaxy contribution in the SN ligh t curve , we em ploy a modified 
version of the ISIS image subtraction routine developed by lAlardl (120001 ) . We revisited 
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the XRF 060218 field nine months after our original observations to obtain deep template 
imaging of the host galaxy when the SN contribution was negligible. In all, we obtained an 
effective exposure of roughly 5.5 hours in the J, H, and K s bands after co-adding several 
nights of imaging of the host galaxy. Using this late-time template, we were able to produce 
residual images containing the measured flux contribution minus the host galaxy, leaving 
only the SN flux. In order to quantify the error involved in this subtraction process, we 
inserted a series of fake stars into the original mosaics and tracked the changes in their 
known magnitude through each subtraction. The resulting variations in the fake star light 
curves are then accounted for in the final host subtracted SN light curve. We find that the 
host galaxy's flux contribution to the uncorrected SN light curve in most of the subtractions 
is roughly equal to the host magnitude found in the late time template imaging, indicating 
that the host galaxy can be treated as a point source at the resolution of our images. This 
allows us to subtract the host contribution in catalog space in order to reduce the scatter in 
the final SN light curve. 



3. Results 

Figure 2 shows our host-subtracted JHK S light curves for SN 2006aj over the course of 
the first 16 days since the burst plotted in the observer frame, with the H and K s bands shifted 
by 0.9 and 1.5 mag respectively for clarity. Each data point reflects the median average of 
all observations taken in a single night of observations. The error bars represent the 1 
cr standard deviation about this median magnitude added in quadrature to the individual 
photometric measurement errors of a single mosaic. On nights in which only one mosaic was 
available, only the photometric error was used, which is most likely an underestimate of the 
true (statistical + systematic) error. 

The characteristic broad peak of SN Ic events can be seen in all three bands, where the 
light curves flatten between day 1 and 13, w i th the K s band peak being more uncertain 



due to a gap in the available data. I Cobb et al.l (120061 ) has reported on J band data peaking 



around day 17, just beyond our observability range, so we may miss the peak magnitude in 
J by a day. To estimate the peak magnitudes of our data set, we employ a spline fit to the 
light curve of each individual filter and measure apparent peak magnitudes of J: 16.76 ± 
0.03, H: 16.65 ± 0.05, and K s : 16.39 ± 0.21. Our estimate of th e J band peak is consistent 



with the peak J band measurements made by ICobb et al.l (120061 ) of the the SN and host of 



16.65 ± 0.06 mag at day 17, one day after our last observation. Both the if and K s bands are 
expected to peak after the J band, so the above values can only act as upper limits to the HK S 
peak magnitudes. All three PAIRITEL bands appear to transition to the plateau phase at 
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roughly the same time, with the J band possibly preceding the H and K s bands. The Galactic 
l ine-of-sight extincti on values for each of the bands, Aj = 0.112, Ah = 0.071, A^ a = 0.045 
( ISchlegel et al.lll998l ). have been subtracte d from the apparent pe ak magnitudes. Extinction 
from the host galaxy has been estimated by lGuenther et al.l (120061 ) through the measurements 
of the equivalent width (EW) of the Na I D lines along the line of sight. They find that most 
of the extinction is attributable to our Galaxy with a visual Galactic extinction of Ay = 0.39 
± 0.02 mag and a corresponding visual host frame extinction of Ay = 0.13 ± 0.01 mag. If 
we assume a ratio of total-to-selective extinction of Ry = 3.1 (ICardelli. Clayton, fc Mathis 



19891 ). we estimate near infrared extinction values from the host galaxy of roughly Aj = 
0.035., Ah = 0.022, Ak s = 0.015 mag, which is considerably less than the e xtinction due to 



the l i ne of sight through our own Galaxy. Because the well- studied hosts (|Vreeswiik et al. 



2004; Jakobsson et al 



20031 : ISavaglio fc Fall! 12004 ; IWatson et al.ll2006l : Butler et all 120061 ) of 



both nearby and cosmological GRBs appear to have a lower dust-to-metal ratio than the 
Milky Way and probably flatter extinction laws, we also consider total-to-selective extinction 
values of R v = 2.0 (5.0), which yield roughly Aj = 0.035(0.035), A H = 0.024(0.022), 
Ak s = 0.015(0.015) mag. None of these values deviates significantly from the Ry = 3.1 
assumption and the effects on the final light curves are within the error of our photometry. 
At a distance of ~ 141 Mpc (z = 0.0335, H = 72 km s -1 , Vt m = 0.3, fi A = 0.7) the peak 
absolute magnitude in J and the upper limit in if and K s , taking into account both sets of 
extinction values, comes to J: -19.02 ± 0.03, H: -19.13 ± 0.05, and K s : -19.39 ± 0.24. In 
this case, these values are not k-corrected due to the low redshift of the event. Photometry 
from the deep late time images gives the host galaxy's JHK S apparent magnitudes J: 18.99 
± 0.16, H: 18.52 ± 0.22, and K s : 18.73 ± 0.34, uncorrected for Galactic extinction. This 
represents 33.6%, 36.7%, and 27.8% of the J, H, and K s flux contribution to the SN light 
curve at early times, respectively. 



It should be noted that fits to the early X-ray data (jButlerll2007l ) collected by the Swift 
spacecraft show evidence for ex cess X-ray absorption bey ond that due to our Galaxy alone 

Dickev fc Lockmanlll99oh: Nn = 2.9 ± 0.5 x 10 21 cm" 2 , at 



A^.Galactic = 1.11 X 10 21 CHI" 2 

z = 0.033. Galactic extinction ( 



Cardelli. Clavton. fc Mathislll989h and the Galactic A V -N H 
relation ( Predehl fc Schmitt 19951 ) would then imply the following absorption magnitudes 
in addition to the absorption by the Galaxy: Aj= 0.51 ± 0.07, A H = 0.33 ± 0.04, A K = 
0.21 ± 0.03. These extinc tion values disa g ree w ith the values calculated through the use of 
the EW of the Na I D by iGuenther et al.l (120061 ) and highlights the uncertainty in applying 
the host galaxy extinction correction. Thus the X-ray inferred results are not subtracted 
from our photometry stated above but should be considered as upper limits to the possible 
extinction from the host galaxy. We briefly discuss possible explanations to these extinction 
disagreements in §3. 
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We compare the observer frame NIR light curve to the light curve behavior at optical 
wavelengths in Figure 3. Here we plo t our JHK S photometry along with UBVRI measure- 
ments made by IMirabal et al.l (l2006bl ) for the first 26 days of the event using the 1.3m and 
2.4m MDM telescope s. Each of the i ndividu al light curves is connected with a cubic spline 
for clarity. The final IMirabal et al.l (l2006bl ) photometry that we adopt for our light curve 
comparison is host-subtracted and extinction-corrected assuming apparent host magnitudes 
of U = 20.10, B = 20.41, V = 20.09, R = 19.91, I = 19.54 and galactic extinction of 
A v = 0.77, A B = 0.61, A v = 0.47, A R = 0.38, and Ai = 0.28. The previously observed 
trend of broader light curves at longer wavelengths can be seen when comparing the turn 
over time in the shorter wavelengths compared to our JHK S measurements. The detailed 
modeling of the light curve behavior and the time-dependent radiative transfer calculations 
req uired for a fu ll understanding of this trend have been accomplished for the case of SNe 
la (IKasenl 120061 ). but have yet to be done for SNe Ib/c and are beyond the scope of this 
paper. This increase in the NIR flux contribution at late times can be nicely seen in the 
color diagram shown in Figure 4. Here we display our host subtracte d and e xtinction cor- 
rected J band photometry minus the UBVRI daXa. from IMirabal et al.l (j2006bl ) and our HK S 
PAIRITEL measurements. The gradual steepening of the color differences is a function of 
AA e ff between the two filters and quantifies the relative difference between the light curve 
decay rates at shorter wavelengths with respect to the NIR. When fit with a linear function 
such that Am = Amo — st, the individual color difference slopes s in Figure 4 are: J—K s : 
-0.020 ± 0.012, J-H: -0.011 ±0.008, J- 1: -0.014 ± 0.006, J-R. -0.032 ± 0.005, J-V: 
-0.048 ± 0.005, J-B: -0.097 ± 0.0046, J-U: -0.143 ± 0.007 mags/day. 

The color evolution properties can also be seen in the spectral energy distribution (SED) 
shown in Figure 5 constructed using this same broad band UBVRIJHK S photometry . We 
converted the UBYRUHK ^ photometry to mono chromati c flux values u s ing Joh nson- Morgan 
/^efr flux zeropoints from iFukugita et al.l (119951) for th e IMirabal et al.l (j2006bl ) data set and 



the 2MASS f ve $ flux zeropoints given in ICohen et all J2003h for our PAIRITEL data. The 
figure shows several SED curves ranging from roughly 1.5 to 12 days after the Swift trigger as 
measured in the observer frame of the host galaxy; the individual SEDs are connected with a 
cubic spline for clarity. From the plot it can be seen that the shorter wavelength contribution 
to the bolometric flux of the SN emission is much higher near the onset of the event and 
steadily decreases after day 7 in the host frame (day ~ 9 in the observer frame). As a result, 
a broadening of the SED profile can also be seen at late times as the NIR contribution to 
the bolometric flux increases, consistent with the conclusions drawn from Figure 4. A much 
more detailed analysis and discussion of the evolution of the XRF 060218/SN 2006aj SED 
using our PAIRITEL observations in conjunction with optical UBVr'i' photometry taken 
with the Mt Hopkins 48-inch will be covered in Paper II (Modjaz et al., in prep). 
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It is difficult to quantify the amount, if any, of the infrared emission over the first 16 
days that is a result of light echos, a scenario in which the dust in the vicinity of the burst 
progenitor absorbs and then reradiates the optical and UV emission from the explosion. This 
would occur at a distance greater than the radius R c at which the dust grains are destroyed 
due to sublimatio n, which when assuming a peak luminosity which is 70% of 1998bw gives a 
rough R c ~ 8-9 pc lWaxman fc Draind (120001 ) . If the dust at this distance has a characteristic 
equilibrium temperature of ~2300 K (the temperature above which the grains are destroyed), 
then the resulting reradiated light would peak at v ~ 2(1 + z)£tm, where z is the redshift of 
the event, which at low z corresponds roughly to the K band IReichartl (120011 ). Although the 
prompt UV and optical emission along the line of sight could contribute infrared emission 
on any time scale, the majority of the reradiated light is expected to be delayed as emission 
arrives from higher latitudes. We conclude that this effect, if present, contributes very little 
to our overall photometry primarily because the large sublimation radius would make the 
peak in the reradiated light occur beyond our last observation on day 16. Furthermore, a 
time varying excess in the K s would be expected if this effect were significant, which is not 
reflected in the SED. 



4. Discussion 



There is no significant evidence that any of the individual NIR JHK S light curves ob- 
served with PAIRITEL peaks significantly faster than the others, although the H and K s 
light curves are expected to peak after the J band which would occur outside of our obser- 
vation window. Furthermore, the rise profile in the individual J, H and K s bands before 
the peak appears to be largely independent of wavelength, with the light curve in each filter 
exhibiting roughly the same slope. This would make the chromatic rise time pr o pertie s of 
SN 2006aj consistent with bolometric rise time results reported by lYoshii et al.l (120031 ) for 
SN 2002ap and points to little or no spectral evolution in the NIR in the days preceding 
peak brightness. This can be seen in the color difference plot shown in Figure 5, which 
shows no significant trend in color evolution between the infrared wavelengt hs during the 



rise portion of the burst. Furthermore, the sharp Fband peak reported by (IModjaz et al. 



2006 



Sollerman et al.ll2006l ; iMirabal et al.ll2006al ) points to a relatively large increase in the 
NIR contribution to the total flux output after day ~ 10. 

The disagreements between the X-ray and optically inferred extinctions are striking but 



not uncommon in the literature (e.g., Campana et all2006l ; iGuenther et al.ll2006l ; ISollerman et al. 



20051 ; IWiersema et al.l 120071 ; IWatson et al.l 120071 ). As mentioned above, GRB host galaxies 
tend to have lower dust to metals ratio than, and probably extinction laws that deviate from 
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that of the Galaxy, possibly explainin g the discrepancy bet ween the Na I D and Nh deter- 
mined extinction values. Furthermore, IWatson et al.l (120071 ) has reported on a disagreement 
between the Nh column densities as inferred from Lya absorption to the metal column den- 
sities from soft X-ray absorption, which they speculate is because the two diagnostics do not 
probe the same gas environments. Their results would imply that the Nh column measured 
from the early X-ray data would have the effect of overestimating the extinction that would 
be inferred through the use of the Ay-Nn relation, possibly explaining the disagreement 
described in S 3. Therefore the exact host extinction could be the effect of a complex picture 
which is still not completely understood. Fortunately, this disagreement on the host galaxy 
extinction affects our NIR data far less than those at shorter wavelengths highlighting the 
importance of NIR observations. 

If we consider day 13 as being the onset of the NIR plateau phase of SN 2006aj, then 
the broad J band peak would occur severa l days after the pea k in the V band, which is 
estimated at 9.7 days in the rest frame by iModjaz et al.l (120061 ) . This would make our J 
band rise time consistent with the previously observed trend in most SNe, in which the 
optical peak precedes the infrared pe ak, as was the case w ith G RB 980425 /SN 1998bw and 
SN 20031w, both of which peaked 1.6 JGalama et alill998h and 5 JMalesani et alil2004h days 
later in / band than in V band respectively. The observed flattening of the NIR as early as 
day 10 makes SN 2006aj the fastest of the XRF/GRB-SNe observed to date, being substan- 
tially shorter in rise time than both GRB 980425/SN 1998bw and GRB 030329/SN 2003dh 
( Galama et al.fl998 ; Matheson et al. The fast evolution and plateau seen in SN 2006aj 

is reminiscent of the light curve properties of broad-lined SNe Ic (SNe Ic BL), highly ener- 
getic core collapse SNe, that do not show likely associations with GRBs. Two such events, 
SN 2002ap and SN 1997ef, have broad bolometric light curves where the latter is es t imate d 
to peak in ~ 10—12 days, assuming the explosion date inferred by iMazzali et al.l (120021 ). 
Similiarly, o ur peak absolute magnitude of —19 .02 ± 0.03 measured in J and —18.7 Mea- 



sured in V (IModjaz et al.l 120061 ; ISollerman et al.l 120061 ; iMirabal et al.l l2006al ) would make 
SN 2006aj the faintest of the XRF/GRB associated SNe, yet still brighter than the av- 
erage broad-lined SN Ic. This can be seen in Figure 6, where the J and R band light 
curves of XRF 060218/SN 2006aj are plotted along the J and R band light curves of 
GRB 980425/SN 1998bw and SN 2002ap. 

Both the peak luminosity and rise time properties of SNe light curves are generally gov- 
erned by the amount of 56 Ni produced in the explosion, leading many authors to conclude 
that the fast evolving and underluminous nature of SN 2 006aj points t o smalle r amounts 
of 56 Ni compared to other XRF/GRB-SNe. Modeling bv IMazzali et all J2006ah estimates 
the 56 Ni mass produced by SN 2006aj to be roughly O.2M , compared to SN 1998bw and 
SN 2003dh which are thought to have produced roughly 0.38— 0.45M and 0.45— 0.65M re- 
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spectively (IMazzali et al.ll2006bl ). Yet, this low estimate of ^ 6 Ni for SN 2006a i is still higher 
than the ~ O.O9M estimated for the SN 2002ap SN Ic BL jFolev et al.ll2003h which reflects 
the amount typical produced by normal core-collapse SNe such as SN 1987A and SN 19941. 

SN 2006aj was also peculiar in several other aspects. First, the measured peak spectral 
energy of the gamma-ray emission E p k = 4.9 keV, where E p k is the max of the vF v spectra 
and hence where most of the gamma-ray energy is radiated, is extremely soft. This is in 
comparison to the GRB counterparts of SN 1998bw, 2003dh and 20031w, which all had 
E p k values of roughly 55 , 79, and 159 keV respectively. Furthermore, the modeling done 
by IMazzali et al.l (j2006al ) estimates an explosion kinetic energy of Ek ~ 2 x 10 51 erg and 



total ejected mass of M e j « 2M Q , both of which are lower than the values estimated for 
the other SNe with associated GRBs, but yet on the high end of the distribution of these 
values when compared to normal SN Ic events. Considering these relatively low E p k, Ek, 
and M e j values along with the low luminosity and fast rise time, it becomes clear that 
GRB 060218/SN 2006aj was indeed a peculiar type of event. It is quite likely that SN 2006aj, 
with its low Ek and M e j values, is an intermediate type of event. Extreme in many ways 
when compared to broad-line core-collapse SNe like SN 2002ap and 1997ef which emit little or 
no high e nergy emission, but rather weak when compared to previously observed XRF/GRB- 
SNe. |Lj (120061 ) has recently quantified a previously observed trend that correlates a GRB's 
E p k value and the peak luminosity of its SNe emission, and hence the produced 56 Ni mass. 
Using this correlation, he estimates that the E p k values for SN 2002ap and 1997ef, if they had 
possessed an associated GRB, would have occurred in the UV, far below the range of what is 
considered a GRB or XRF. This of course is based on the assuming this correlation is real and 
not an artifact of observational biases or source evolution. All of the SN 2006aj properties 
discussed above place it between the previously observed GRB-SNe sample and these low- 
luminosity and low estimated E pk broad-lined events, possibly pointing to a continuous 
distribution of SN events that emit at increasingly higher energies. 

One possible explanation for such a continuous distribution of gamma ray energies 
is that the soft-spectrum and low-luminosity ev ents like XRF 060218 are due to off-axis 
observations of typical long-duration GRB-SNe (jYamazaki et al.1 120031 ). Although this off- 
axis scenario was originally involked to explain the peculiar high e nergy proper ti es of GRB- 
SNe like SN 1998bw, recent numerical simulations performed by iMaeda et al.l (120061 ) have 
shown that the viewing angle of asymmetric explosions would also have an effect on the 
observed luminosity and time of peak of the associated SNe. Through the use of 3D Monte 
Carlo simulations, the authors find that asymmetric SNe would appear to peak earlier when 
viewed at small 9 from the z-axis of the explosion, this being due primarily to low and 
extended 5 6Ni densities in that direction. The resulting luminosity as seen in the z-direction 
is then boosted correspondingly because the photons can only diffuse out in this direction. 
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This small 9 scenario would then predict fast peaking, high luminosity, and high E p k events 
which is inconsistent with the parameter distribution seen in XRF 060218 which has a 
distinctly fast time to peak but a remarkably low luminosity and E p ^. Another difficulty 
of the off-axis interpretation is that it would require a very large intrinsic total energy of 
roughly ~ 10 53 ergs, which would be incon sistent with the relati vely low kinetic energy 
measurements made from radio observations (jSoderberg et al.ll2006l ) at late times when the 
jet collimation and/or SNe asymmetry should no longer be significant. Furthermore, the 
probability of seeing a jet off-axis goes down significantly with increa sing 0j resu l ting i n 
inconsistent rate ratios between high and low energy events as noted by lCobb et al.l (120061 ). 
Therefore, although off-axis models can be made to successfully explain the observed range 
of GRB/XRF energetics, the range of SN properties such as Ek, M e j, the measurements 
of which are generally independent of collimation, suggest that viewing angle effects alone 
cannot account for the full diversity of XRF/GRB-SNe properties. 

On the other hand, the distribution of the intrinsic properties of the progenitor star 
that would lead to such a continuous distribution in the observed XRF/GRB-SNe properties 
is also not entirely clear. A qualitative correlation between the estimated progenitor mass 
and t he resulting explosiv e energy seen among several broad-line and normal core-collapse 
SNe (IMazzali et al.l |2002| ) points to a high initial stellar mass along with the progenitor's 
rotation rate, both of which would govern whether the jet associated with the GRB breaks 
out of the progenitor, as likely candidates. Nonetheless, future broad-band observations of 
intermediate events like SN 2006aj and SN 2002ap may be crucial to the examination of this 
distribution. 



Finally, we note that the iGorosabel et aJj (120061 ) have reported evidence of a possible 
circularization of the explosi on geometry, which l ike previous XRF/ GRB-SN is thought to 
be aspherical at early times (IMazzali et al.ll200ll ; iHjorth et al.ll2006l ). The authors show a 
rotation in the polarization angle of ~90 degrees at 14 < t < 39 days whic h they interpreted 
as a circularization of the previous highly aspherical geometry. Similarly, iPian et al.l (120061 ) 
reported a change in the decay slope of the photospheric expansion velocities vs. time 
as determined through modeling of the spectra at the various epochs, which could also 
be interpreted as evidence for circularization. This circularization time of t ~ 14 roughly 
corresponds with the onset of the plateau phase in our NIR observations of SN 2006aj. One 
could speculate that the fast rise observed in SN 2006aj is a result of fast and aspherical 
outflow which becomes spherical after t ~ 14 days, creating the observed plateau phase. 
Although, the fact that a general broadening of SN light curves in the NIR are seen in other 
types of (presumably spherical) SN events would point to a mechanism common to all SNe, 
such as radiative transfer and/or the ionization changes in the SN ejecta (lKasenll2006l ) . rather 
than the geometry of the explosion. Future NIR observations on the existence, or lack thereof, 



of a NIR plateau in normal (non-broadlined) SN Ic along with late-time spectral modeling 
and polarization measurements should shed light on any potential correlation between the 
expansion geometry an the resulting NIR light curve. 
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Figure Captions 



Fig. 1. - Color composite finding charts of the host galaxy (seen in SDSS pre-imaging; 
left) and SN 2006aj (seen in PAIRITEL; right). The SDSS image is made using r, i, z-band 
images (Cool et al. 2006) and the IR image using J, H, K s from stacked mosaics obtained 
on 3 March 2006 UT, near in time to the peak of the optical emission of the SN. North is 
up and East is to the left. 

Fig. 2 - The J band (blue circles), H band (green circles), and Ks band (red circles) 
aperture photometry light curve of XRF 060218/SN 2006aj after subtraction of the host 
contribution. The if and K s band light curves have been shifted by —0.9 and —1.5 magni- 
tudes respectively for clarity. The three light curves have been corrected for Galactic and 
host line of sight extinction. The error bars represent the quadrature sum of the photometric 
errors associated with all the mosaics produced in a single night of observations plus the 1 
a standard deviation about their median magnitude. A flattening of the light curves in all 
wavebands can be seen by Day 16. 

Fig. 3 - Our near infrared PAIRITEL light curves plotted along with UBVRI data 
taken from IMirabal et al.l (l2006bl ). A steady progression of wider light curves at longer 
wavelengths is clearly seen. 

Fig. 4 - Color diagram showing the subtracti on of the PAIRITEL J band and H and 
K s band data as well as the IMirabal et al.l (l2006bl ) UBVRI photometry. The slope of the 
color differences becomes steeper with increasing AA e g, showing a large increase in the NIR 
contribution to the total bolometric flux output of the SN at later times. 

Fig. 5 - A br oad spectral energy di stribution diagram produced by converting the JHK S 
band data and the IMirabal et al.l (l2006bl ) UBVRI photometry into monochromatic flux values 
and plotted vs the effective frequency of that band pass, corrected for the redshift of the 
host galaxy. The NIR contribution to the total energy represented by the SED can be seen 
to increase at later times. 



This preprint was prepared with the AAS IATgX macros v5.2. 
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Fig. 6 - A comparison of the R and J band light curve pro perties of XRF 060 218/SN 2006 aj 



(PAIRITEL), GRB 980 425/SN 1998bw jGalama et al.lll998h . and SN 2002ap JFolev et al. 



20031 : lYoshii et al.1 l2003h . We us e the explosion date of January 28.9 UT (JD = 2452303.4) 
inferred by iMazzali et all hoO$ \ as the t = time for SN 2002ap. 
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2006aj Multicolor Light Curves 
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2006aj Color Differences 
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2006aj Spectral Energy Distribution 
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Table I. Host Subtracted Photometry of XRF 060218/SN 2006aj 



Days 


J band 


if band 


K s 


band 


Seeing 


Air Mass 


since GRB 1 


mag 2 


mag 


2 


mag 


2 


(arcsec) 


(sec z) 


1.99 


18.31 


± 0.09 


17.97 


± 


0.21 


17.72 


± 


0.48 


2.93 


1.10 


2.98 


18.18 


± 0.09 


18.04 


± 


0.33 


17.26 


± 


0.15 


2.64 


1.11 


4.01 


17.85 


± 0.19 


17.89 


± 


0.17 


17.10 


± 


0.10 


2.76 


1.20 


4.98 


17.58 


± 0.02 


17.59 


± 


0.13 


16.91 


± 


0.20 


2.77 


1.22 


5.98 


17.25 


± 0.05 


17.34 


± 


0.13 


16.63 


± 


0.42 


2.89 


1.32 


6.96 


17.15 


± 0.07 


17.15 


± 


0.11 


16.69 


± 


0.10 


2.85 


1.32 


7.97 


17.12 


± 0.04 


17.13 


± 


0.30 


16.28 


± 


0.11 


2.92 


1.43 


8.96 


16.96 


± 0.06 


16.89 


± 


0.11 


16.35 


± 


0.16 


2.95 


1.29 


9.98 


17.14 


± 0.07 


16.79 


± 


0.04 








2.85 


1.29 


12.98 


16.88 


± 0.04 


16.85 


± 


0.03 


16.32 


± 


0.21 


1.45 


1.21 


13.99 


16.97 


± 0.02 


16.73 


± 


0.05 








2.82 


1.22 


14.96 


16.89 


± 0.03 


16.83 


± 


0.05 


16.48 


± 


0.09 


2.64 


1.39 


16.00 


16.93 


± 0.05 














2.83 


1.37 



1 Observer frame 
Uncorrected for extinction 
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Table 2. XRF 060218/SN 2006aj Properties 



Property 


J band 


H band 


K s band 




mag 1 


mag 1 


mag 1 


Peak mag 


16.76 ± 0.03 


< 16.65 ± 0.05 


< 16.39 ± 0.21 


Peak Mag 2 


-19.02 ± 0.03 


< -19.13 ± 0.05 


< -19.39 ± 0.24 


Host mag 


18.99 ± 0.16 


18.52 ± 0.22 


18.69 ± 0.34 


Agalactic 


0.112 ± 0.003 


0.071 ± 0.002 


0.045 ± 0.001 


A 3 


0.035 ± 0.003 


0.022 ± 0.002 


0.015 ± 0.001 



^ost subtracted and extinction corrected values 
2 z = 0.0335, H = 72 km s" 1 , tt m = 0.3, fi A = 0.7 
3 R V = 3.1 



Table 3. XRF 060218/SN 2006aj Absolute Magnitudes 



Days 


J band 


H band 


K s band 


since GRB 1 


9 3 

mag ' 


9 ^ 

mag ' 


9 ^ 

mag ' 



1. 


92 


-17. 


.64 


± 


0. 


.09 


-17. 


.92 


± 


0. 


.21 


-18 


,13 


± 


0.49 


2, 


.88 


-17 


,76 


± 





.09 


-17. 


.86 


± 


0. 


33 


-18. 


.25 


± 


0.15 


3. 


.88 


-18. 


.08 


± 





.19 


-18. 


00 


± 


0, 


.17 


-18. 


.45 


± 


0.10 


4. 


.81 


-18. 


.34 


± 





.03 


-18. 


.29 


± 


0, 


.13 


-18 


.69 


± 


0.22 


5. 


.79 


-18. 


.66 


± 





.06 


-18. 


53 


± 


0, 


.14 


-19. 


.02 


± 


0.48 


6. 


73 


-18. 


.76 


± 





.08 


-18. 


.72 


± 


0, 


.13 


-18. 


.96 


± 


0.11 


7. 


.71 


-18. 


.79 


± 





.05 


-18. 


.74 


± 


0. 


33 


-19. 


.42 


± 


0.14 


8. 


.67 


-18. 


.95 


± 





.07 


-18. 


.98 


± 


0. 


.12 


-19 


,35 


± 


0.19 


9. 


.65 


-18. 


,77 


± 





.07 


-19. 


.08 


± 


0. 


.04 










12. 


.56 


-19. 


.02 


± 





.05 


-19. 


.01 


± 


0, 


03 


-19 


.38 


± 


0.24 


13. 


.54 


-18. 


.94 


± 





.02 


-19. 


.13 


± 


0. 


.05 










14. 


.47 


-19. 


.01 


± 


0. 


.03 


-19. 


04 


± 


0. 


06 


-19. 


,19 


± 


0.11 



iRest frame of GRB 

2 Host subtracted and extinction corrected values 
3 z = 0.0335, H = 72 km s" 1 , Vt m = 0.3, fi A = 0.7 



